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INTERPLANETARY MEDIUM DATA BOOK - SUPPLEMENT 3 
1977 - 1985 


ABSTRACT 

The updating of the hourly resolution near-Earth solar wind data 
compilation is discussed. Data plots and listings are then presented. 

In the text, the time shifting of ISEE 3 magnetic field and plasma data, 
using corotation delay is explained in detail. Normalizations of 
IMP (MIT) , ISEE 3, and ISEE 1 temperatures and densities to equivalent 
IMP(LANL) values are also discussed. Hie levels of arbitrariness in 
combining data sets, and of random differences between data sets, are 
elucidated. 


vii. 




Introduction 


In previous issues of this series, hourly averaged interplanetary 
magnetic field and plasma parameters have been listed and plotted. These 
parameters have come from a number of spacecraft in the near-Earth solar wind. 
Data for 1963-1975 were contained in the Interplanetary Medium Data Book 
(NSSDC/WDC-A-R&S 77-04, 1977), data for 1975-1978 were published in the 
Interplanetary Medium Data Book - Supplement 1 (nssdc/wdc-a-r&S 79-08, 1979) 
and data for 1978-1982 were published in the Interplanetary Medium Data Book 
- Supplement 2 (NSSDC/WDC-A-R&S 83-01, 1983). This third supplement 
represents an extension of the earlier compilations to 1985. Supplement 3 
supersedes Supplement 2 in terms of both the data and the discussion. It 
begins in 1977 in order to incorporate ISEE 1 data, and corrects previously 
announced errors published in Supplement 2 which affected values of B y (GSM) 
and B z (GSM) in 1980 to 1982 and some By (GSE) values in raid-1980. Because of 
its size, this supplement has been published in two volumes. The first volume 
(Supplement 3) contains descriptive information and 27-day plots of various 
parameters. Listings of selected parameters are contained in Supplement 3A. 
Copies of all the books in this series, as well as the parent OMNI tape from 
which the listings and plots are generated, are available from the National 
Space Science Data Center (NSSDC). These data (from 1976 onward) are also 
accessible online through the NSSDC VAX. Access procedures are described 
subsequently . 


Data Sources for this Supplement 

Data for this supplement come from the IMP 8, ISEE 3 and ISEE 1 
spacecraft. IMP 8 was launched on October 26, 1973, and is in a low 
eccentricity, ~30Re x 40 Re, 12.5-day geocentric orbit. IMP 8 is in the solar 
wind for a duration of 6 to 8 days per orbit. ISEE 3 was launched on August 
12, 1978. Until mid- 1982 it orbited the Earth-sun libration point 
approximately 240 Earth radii sunward of the earth, ranging up to 100 Re from 
the Earth-sun line. ISEE 3 was then redirected into the Earth's geotail. It 
sampled the solar wind for a large portion of 1982 and for portions of 1983. 
Launched on October 22, 1977, ISEE 1 has a highly elliptical orbit with an 
apogee of 23 Re and a period of about 57 hours. Usually ISEE 1 only enters 
the solar wind from August through December where it can spend up to 30 or 
more hours of each orbit in the interplanetary medium. As of early 1986, IMP 
8 and ISEE 1 continue to provide near-Earth data. 

The six data sets folded into this compilation are identified in Table 1. 
Some of the plasma data attributed to Los Alamos IMP 8 actually were provided 
as a mix of Los Alamos IMP 6, IMP 7 and IMP 8 data. All ISEE 1 and IMP 8 
plasma parameters are based on ion measurements. ISEE 3 plasma parameters are 
based on ion measurements through day 48, 1980 (at which time the ion portion 
of the instrument failed), and electron measurements after this time. When 
both ISEE 3 and IMP 8 (field or plasma) data were available for a given hour, 
the IMP data were used in this compilation due to IMP'S greater proximity to 
the Earth. When both MIT and LANL IMP 8 plasma data were available, the MIT 
data were chosen for historical reasons. ISEE 3 plasma data were chosen 
preferentially to ISEE 1 plasma data because of the availability of additional 
statistical and flow direction parameters. 
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TABLE 1. Data Sources for this Supplement 


Spacecraft 


Principal Investigator Time Span (YR/DOY) 


Magnetic Fields 
Plasma 


IMP 8 N. 
ISEE 3 E. 
IMP 8 H. 
IMP 8 S. 
ISEE 3 S. 
ISEE 1 S. 


F. Ness (GSFC) 
J. Smith ( JPL ) 
S. Bridge (MIT) 
J. Bame (LANL) 
J. Bame (LANL) 
J. Bame (LANL ) 


73/302 85/091 

78/225 83/365 

73/334 85/107 

71/076 84/357 

78/228 82/279 

77/303 79/365 


* Key scientists associated with the reduction of these data include 

R. P. Lepping, B. T. Tsurutani, J. D. Sullivan/ A. J, Lazarus, J. A. Gosling, 

W. C. Feldman, and R. D. Zwickl. 


Systematic and Random Differences Between Data Sets 

The principal task in interspersing data from different sources is to 
make the data sets as compatible as possible. Thus, a data compilation is 
desired in which small amplitude, long term changes, and larger amplitude, 
short term changes may be studied with relatively high confidence levels. 
Absolute values of parameters, and small amplitude, hourly changes are less 
accurately determined, especially when more than one data source is involved. 

There may be random differences and systematic differences between data 
sets. Systematic differences result from calibration errors in one or both 
data sets* In this context, "calibration" implies the whole process of 
forming the electrical output of a sensor or sensors measuring some 
physical parameter to a declaration of what the physical parameter value is. 


Random differences occur at least in part because two instruments may be 
measuring at different times or places, in the presence of temporal evolution 
or spatial gradients in the parameter being measured. In the case of hourly 
averages, differing parts of an hour may have been measured ky two 
experiments. Random differences may also arise from the use of the differing 
instrumentation and data analysis technique, which has been discussed by 
Neugebauer ( Sp . Soi. Rev, , 33* 127, 1982). 

When a regression analysis is performed between two data sets of 
"simultaneously measured" parameter values, the deviation of the regression 
line from Y = X gives a measure of the systematic difference between data 
sets, and the scatter of data points about the regression line gives a me as ure 
of the random difference. For the case of ISEE 3, data must first be time 
shifted to expected Earth arrival times. If this is not done, the random 
differences found by regressing IMP 8 and ISEE 3 parameters are not 
irreducible. This is because the ISEE 3 to IMP 8 transit time is comparable 
to the hourly time resolution used in this conpilation. 

For those parameters for which systematic differences are comparable to 
or greater than the random differences, a cross-normalization of data sets is 
performed. It has been found that only density and temperature need to be 
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cross-normalized. That is, IMF parameters and flow speed had systematic 
differences significantly less than random differences. The decisions as to 
which data set to normalize to which other data set involved a certain 
arbitrariness, which was discussed in the original Data Book and Supplement 1. 

Time Shifting of ISEE 3 Data 


Either or both of two approaches to time-shifting ISEE 3 fine scale data 
may be expected to yield statistically irreducible random errors. (There are 
random differences between two sets of observations and random errors between 
either set of observations and the "true" values.) These approaches are 
corotation and convection. Consider initially the first of these. The 
corotation delay equation is: 


r 



In this equation, t is the time shift, X and Y are the ecliptic plane 
projections of the geocentric ISEE 3 position vector along and across the 
Earth-sun line, V is the measured solar wind speed (assumed radial), Rfl is the 
equivalent speed of solar rotation at 1 AU (-428 km/s), and Ve is the orbital 
speed of the Earth (~30km/s). 

Use of the corotation delay equation would be completely correct if the 
constant-phase surfaces of all solar wind variations were normal to the 
ecliptic plane and aligned with the ideal spiral interplanetary magnetic field 
(IMF) direction, and if solar wind variations had no temporal evolution on the 
ISEE-to-Earth transit time scale nor spatial gradients transverse to the 
Earth— sun line on scales less than —100 Re* None of these required 
conditions is satisfied by all solar wind variations. Thus, random errors are 
an inevitable part of the inference of near -Earth solar wind parameter values 
from ISEE 3 data and the corotation delay equation. 

A study of these errors was begun by comparing "simultaneous" 
5-min-averaged ISEE 3 and IMP 8 field and plasma data. Simultaneity here 
implies that the IMP 8 time is within 2. 5 min of the corotation-shifted ISEE 3 
time. In this analysis, the X and Y of the corotation delay equation relate 
to the ISEE-IMP separation vector rather than to the geocentric ISEE position 
vector . 

Figure 1 shows a scatter plot of simultaneous, 5-min averaged ISEE 3 and 
IMP 8 measurements of the Z(GSE) component of the IMF. These data were taken 
between day 36 and 66 of 1979, when ISEE 3 was within 50 R e of the Earth -sun 
line. There are 1867 data points included. The regression line determined by 
minimizing the sum of squares of perpendicular distances between data points 
and regression line is shown on the figure. That the regression line is not 
Y “ X relates to systematic differences and will be addressed in the next 
section. Of principal interest here is the scatter of points about the 
regression line. The rms perpendicular distance between data points and 
regression line (labelled RMS PERP DIST OF PTS ABT LINE on the plot and called 
the "scatter function" for convenience in the following) is 1.35 nanoTesl as 
(nT). This means that for a given ISEE 3 5-min observation of B z , IMP 8 would 
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Figure 1. Scatter plot and best fit regression line for 
IMP 8 and ISEE 3 B z 5-min values, Feb. 6 - Mar. 7, 1979 


4 



measure, over the five corresponding corotation -shifted minutes, a value of B z 
in the range given by the ISEE 3 value, plus or minus V 2 * 1.35 nT, with a 
probability of ~68%. (The \f2 factor enters because 1.35 nT is a distance 
normal to the regression line, while for the present purpose we need an 
equivalent distance parallel to the ordinate. Regression line slopes are 
close enough to unity that \[2 is an adequate approximation throughout. ) 

Using 5-min data, regression lines and scatter functions have been 
determined for each of seven interplanetary parameters, for each of eight 
(magnetic field) or thirteen (plasma) time periods, for each of three 
approaches to time-shifting ISEE 3 data. (Time periods after day 48, 1980, 
were treated separately since 5-min resolution plasma data were not available. 
The method used to time shift these later time periods is explained toward the 
end of this section. ) The parameters analyzed are field Cartes ian components 
(B x , By, B z ), field magnitude (B m ) , flow speed (V), proton density (N), and 
temperature (T). Because N and T cure distributed more logarithmically than 
linearly, log N and log T were used in the regression analysis. 

The time periods for which analyses were performed were defined by the 
ISEE 3 orbital phase: Y(ISEE) < -50 Rg? -50 Rg < Y(ISEE) < 50 Rg? 50 Rg < 

Y(ISEE). Most intervals so defined had durations between 30 and 60 days, and 
contained between 1000 and 3000 ISEE/IMP 5-min data points. Note that because 
IMP 8 covers the range -40 Rg < Y < 40 Rg during each 12.5 day orbit, there is 
some overlap of Y( ISEE )-Y( IMP) between |y(ISEE)[ < 50 Rg and |y(ISEE)| > 50 Rg 
periods. 

The three time shifts used for ISEE data were corotation (delay equation 
given above), convection [delay = (X ( ISEE )-X( IMP) )/V] , and no shift. 

Convection is expected to yield smaller random errors than corotation if the 
constant-phase surfaces of interplanetary variations are more nearly normal to 
a heliocentric radius vector than they sure aligned with the ideal spiral IMF 
direction, by a statistically significant amount. The no-shift case was 
included to estimate the level of reduction of the random differences produced 
by time-shifting ISEE data. 

Figure 2 summarizes the scatter function information for all physical 
parameters and time periods, for the case of corotation-shifting of the ISEE 3 
data. Note that scatter functions separately computed for individual 
Cartesian components of the IMF have been averaged. Note also that scatter 
functions computed for |Y(ISEE)| < 50 Rg and |Y(ISEE) > 50 Ro intervals are 
distinguished. Lines have been drawn connecting the Y(ISEE)| < 50 Rg 
intervals. 

Several points may be noted in Figure 2. First, there is variability 
from one interval to the next. For the IMF components , the scatter functions 
are larger for the |y(ISEE)| > 50 Rg intervals them for the |y(ISEE)| < 50 Rg 
intervals. This is less clearly true for the field magnitude, and is not true 
for the plasma parameters. These statements are also true when convective 
shifting of ISEE data is considered. This behavior is consistent with the 
transverse scale of plasma parameter variations being greater than 50 to 100 
R e , and with the scale for field component (or direction) variations being 
comparable to or less than 50 to 100 Rg. Field magnitude variation scales 
seem to be intermediate. 
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Scatter Functions, Corotation Shift 


5min. Analysis 



Figure 2. Scatter functions, corotation shift, 5-min analysis 



Note that the scatter functions for field magnitude are significantly 
smaller than scatter functions for field compone nts . This is consistent with 
a preponderance of constant-field - intensity Alfvenic fluctuations in the 
solar wind. 

As averaged over all time periods, the mean scatter functions for the 
five panels of Figure 2 are 1.52 nT, 0.51 nT, 12.4 km/s, 0.084, and 0.108. 
Recall the above discussion of Figure 1 concerning the relation of these 
values to the inference of near-Earth parameter values from ISEE observations. 
Note that the last two of these five mean scatter functions correspond to 
one-sigma uncertainties in near— Earth density and temperatures values of 30% 
and 40%, given ISEE 3 observations. 

Next we compare these results, based on corotation shifting, with the 
results of convective shifting and no shifting. Figure 3 (left and center 
panels) contains histograms of ratios of scatter functions for these latter 
analyses to those illustrated in Figure 2. Note that |y(ISEE)| < 50 Rq and 
|Y(ISEE) | > 50 Rq cases are again distinguished. 

From the convection-to-corotation histogram, it is apparent that very 
similar scatter functions are obtained with corotation and convection. Very 
seldom do the convection and corotation scatter functions differ ty more than 
10% for a given physical parameter and time period. When all time periods and 
parameters are considered, it is found that the convective scatter function 
exceeds the corotation scatter function in 55% of the cases (39 cut of 71). 
Even for the |y(ISEE)| > 50 Rg cases when the corotation and convection 
analyses might be ejqpected to exhibit the greatest differences (owing to the 
greater delay differences), there is no statistically significant difference 
between the corotation and convection cases. From this it may be concluded 
that the distribution of const ant -phase surfaces of interplanetary variations 
is peaked neither along the ideal IMF spiral direction nor perpendicular to a 
heliocentric radius vector. Our choice of time shifting by corotation rather 
than by convection is seen to be a basically arbitrary choice. 

Given this statistical equivalence of convective and corotational time 
shifting of ISEE 3 data, it may be asked whether these shifts give a 
statistically significant improvement (reduction of random differences, or 
scatter functions) relative to not shifting ISEE data at all. From the 
no-shift to corotation-shift histogram of Figure 3 (center panel), it is clear 
that time shifting does indeed make a difference. The extent of the 
difference depends on both the physical parameter and on |y(ISEE)|. 

Time shifting reduces the random differences by anywhere between 0% and 
40% for plasma parameters, and up to a factor of 2 for field parameters. The 
greater improvement in field parameters results from these parameters * smaller 
scale lengths. That is, a greater portion of the random differences between 
(no-shift) simultaneous ISEE 3 and IMP 8 field parameters results from the 
temporal -spatial displacement of these spacecraft than is the case for plasma 
parameters, and this displacement is at least partially compensated for by the 
time-shif ting. 

Time shifting also reduces the scatter functions by a greater amount for 
the | Y( ISEE) | < 50 Re intervals than for the |y(ISEE)| > 50 Re intervals; this 
effect is more pronounced for the field parameters than for the plasma data. 
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This is an effect of transverse spatial gradients in interplanetary variations 
and the existence of a distribution of constant-phase orientations. These 
become increasingly significant as |y(ISEE)| increases, and they are handled 
only imperfectly with a simple time shift which assumes a uniform 
constant-phase orientation. 

All the foregoing discussion has concerned 5-min resolution analysis. It 
is expected that coarser time resolution would yield smaller random 
differences, since finer scale structures which probably contribute 
significantly to the 5-min-resolution random differences would be averaged 
out. This expectation is realized. 

ISEE 3 hourly averages, constructed from corotation-shifted 5-min 
averages, have been compared to corresponding IMP 8 hourly averages. The 
scatter functions determined in this analysis have been compared to those 
obtained in the 5-min resolution analysis. Ratios of hour-based scatter 
functions to 5-min-based scatter functions are also shown in Figure 3 (right 
panel). The process of taking hourly averages of 5-min averages reduces the 
random difference between data sets by up to a factor of 2. There is no 
major dependence on which physical parameter is involved, nor on |y(ISEE)|. 
There is a weak suggestion that the reduction in scatter function is slightly 
greater for flow speed and for field components than for the other 
parameters. 

In Figure 4 are presented the scatter function values themselves, as 
obtained from the hourly resolution analysis. The format and scaling are the 
same as for Figure 2, although the scales have been translated to account for 
the smaller hourly based scatter functions. There were typically between 300 
and 600 pairs of "corotation-simultaneous " ISEE and IMP hourly averages in 
each time period. 

Again, variability is observed from one interval to the next, although 
for most parameters the amplitude of the variability is reduced relative to 
Figure 2. Again, scatter functions for field components are smaller at 
|Y(ISEE)| < 50 Re than at |y(ISEE)| > 50 Re? this |y(ISEE)| dependence is less 
apparent for field magnitude and is not apparent for plasma parameters. 

As discussed in the next section, a single cross-normalization will be 
applied to each parameter for all time. It is of interest to determine random 
differences between ISEE and IMP for each parameter over the full span of the 
simultaneously available data. Figures 5-11 contain scatter plots, regression 
lines and scatter functions for all seven physical parameters, based on all 
available "corotation-simultaneous" hourly averages. 
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PARAMETER : BX - C5E 

TIME SPAN IS: 78224 TO 79178 

NO. OF PTSs 3400 
SHIFT TYPE: COROTATION 

AVERAGES: IMP 0.24 14.14 

■ ISE 0.20 14.22 

BX-IMP - 0.04 110.02) + 0.98(10.01) * 8X-ISEE 
RMS PERP DIST OF PTS AST LINE: 0.97 


■r. 


-m.oo 


ISEE-3 


BX - GSt 


Figure 5. Scatter plot and best fit regression line for 
IMP 8 and ISEE 3 B x hourly averages, ftng. 78 - June 79 


PARAMETER i BY - CSE 
TIME SPAN ISj 78224 TO 79178 
NO. OF PTS: 3400 

SHIFT TYPE: COROTATION 

AVERAGES: IMP -0.15 14.73 

: ISE -0.1S 14.87 

BY-IMP • -0.00(10.02) ♦ 0.97(10.01) » BY-ISEE 
RMS PERP DIST OF PTS ABT LINE: 1.01 


ISEE-3 


BY - GSE 


Figure 6. Scatter plot and best fit regression line for 
IMP 8 and IS EE 3 B„ hourly averages, ftug. 78 - June 79 



21.00 


PARAMETER * BZ - GSE 
TIME SPAN IS: 78224 TO 79178 

NO. OF PTSi 3400 
SHIFT TYPE: COROTATION 

AVERAGES* IMP -0.23 13.55 
« ISE -0.30 13.47 

BZ-IMP - 0.07(10.02) ♦ 1.03(10.01) h BZ-ISEE 
RMS PERP DIST OF PTS RBT LINE: 1.12 


I SEE-3 


BZ - GSE 


Figure 7. Scatter plot and best fit regression line for 
IMP 8 and ISEE 3 B z hourly averages, Aug. 78 - June 79 




B-MRGNI TUDE 


PARAMETER : B-MRGNITUDE 

TIME SPAN IS: 78224 TO 79178 

NO. OF PTS: 3400 

SHIFT TTPE: COROTATION 

AVERAGES: IMP 7.30 13.08 

: ISE 7.42 13.08 

MHG-IHP ■ -0.13110.01) ♦ 1.00(10.00) « MAG- 1 SEE 
RMS PERP OIST OF PTS ABT LINE: 0.48 


I SEE-3 B-MRGNITUDE 


Figure 8. Scatter plot and best fit regression line for 
IMP 8 and ISEE 3 hourly averaged field magnitudes, 

Aug. 78 - June 79 





PARAMETER : SPEED 

TIME SPAN ISs 78228 TO 80049 
NO. OF PTS: 5000 

SHIFT TYPE: COROTATION 

AVERAGESi IMP 417 181 

t ISE 410 ±81 

V-IMP • 5.97 (±0. 051 ♦ l.OOIJO.OO) x V-ISEE 
RMS PERP DIST OF PTS ABT LINE: 8.44 



335.00 410.00 485.00 560.00 635.00 710.00 

I SEE-3 V (KM/S) 

Figure 9. Scatter plot and best fit regression line for 
IMP 8 and IS EE 3 hourly averaged speeds. Aug. 78 - Feb. 80 
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LOG DENSITY 


PARAMETER : LOG DENSITY 

TIME SPAN ISs 78228 TO 80049 
NO. OF PTSs 5000 
SHIFT TYPE: COROTATION 

AVERAGES: IMP 0.82 10.32 

i ISE 0.79 10.32 

LOG-N-IMP - 0.04(10.00) + 1.00(10.00) * LOG-N-ISEE 
RMS PERP DIST OF PTS ABT LINE: 0.07 






xv 




£ •} 


I SEE-3 LOG DENSITY 


Figure 10. Scatter plot and best fit regression line for 
IMP 8 and ISEE 3 logarithms of hourly averaged densities, 
Aug* 78 - Feb. 80 






PARAMETER : LOG TEMP 

TIME SPAN IS: 78228 TO 80049 

NO. OF PTS: 5000 

SHIFT TYPE: COROTATION 

AVERAGES: IMP 4.95 10.30 

: ISE 4.90 10.29 

LOG-T-IMP ■ -0.08(10.01) ♦ 1.03(10.01) * LOG-T-ISEE 
RMS PERP DIST OF PTS ABT LINE: 0.09 




The scatter function information contained on those figures may be summarized 
by the first part of Table 2 below and the statement: 

If an hourly averaged value of X is obtained for parameter Y from remotely 
measured but corotation-shifted ISEE 3 fine-scale data, then the near -Earth 
IMP 8 would measure a value for parameter Y in the range X + Z with a 1 -sigma 
(~68%) confidence level, where the Z functions are given below. 


TABLE 2. Summary of Hourly Z -Function Information 


Y 

5 -min 
for t 

ion -based speeds used 
calculation (1978-1980) 
Z 

Hourly electron-based speeds 
for x calculation (1979-1983) 
Z 

B x 


1.37 nT 

1.63 

By 


1.43 

1.80 

jr 

B z 


1.58 

1.88 

B m 

V 


0.68 

0.75 


11.9 km/sec 

29.3 km/sec 

log N 


0. 10 

0. 14 

log T 


0. 13 

no ion tempera- 
tures available 


Note that the Z factors differ from the scatter functions ty as 
explained previously. Vfe take these values to represent the irreducible 
random differences arising from the interspersal of ISEE 3 and IMP hour 
averages. (The right-hand column in Table 2 is discussed below. ) 

All of the above discussion of time shifting concerned data that were 
collected prior to day 48, 1980, and used originally in Supplement 2. After 
this time, the ISEE 3 plasma data available were hourly averaged parameters 
based on electron observations. The LANL group has carefully normalized their 
densities and velocities to their ion-based data fcy comparing parameters when 
both ion and electron measurements were taken. Nonetheless, we should bear in 
mind that errors associated with electron-based parameters are somewhat higher 
than those associated with ion— based parameters. 

After day 48, 1980, tim e shifting of ISEE 3 IMF data was effected using 
the hourly resolution ISEE 3 speeds. (For late 1982 and 1983 when ISEE 3 IMF 
data were available but ISEE 3 plasma data were not, hourly resolution IMP 8 
speeds were used.) Note that for portions of 1982 and 1983, ISEE 3 was 
actually tailward of the earth and thus time shift values can be positive or 
negative. Using the time-shifted ISEE 3 IMF hourly values, regression lines 
and scatter functions for ISEE 3 vs. IMP 8 were determined for B x , By, B z and 
B m . The scatter functions are summarized in Table 2. 
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IMF scatter functions calculated when only hourly averaged speeds were 
available have slightly higher values than those calculated when 5 minute 
averaged speeds were used for time shifting. This result is not totally 
unexpected especially when one recalls that the hourly averaged speeds are 
based on electron measurements. It is also important to point out that the 
two separate analyses (5-min ion-based speeds and hourly electron- based speeds 
used for time shifting) were performed on different time intervals. The 
scatter functions do vary somewhat from year to year and this may account for 
much of the differences. 

In the absence of 5-min resolution plasma data, time shifting of the post 
day 48, 1980, plasma data involved a much less straightforward procedure than 
that which was done previously. Each ISEE 3 hourly averaged density and speed 
was shifted to the Earth using the observed ISEE 3 speed and the corotation 
delay equation. Typically an ISEE hour (e.g. , 0200-0300) shifted to a 
non-integral Earth hour (e.g., 0245-0345). Earth-time hourly averages (e.g., 
for 0300-0400) were then built as averages of ISEE values which shifted into 
the desired Earth hour, weighted by their relative contribution to that hour. 

In order to determine the value of the weighted average approach to time 
shifting, a regression analysis was done between IMP 8 plasma data and ISEE 3 
plasma data which were not time shifted. Scatter functions calculated for the 
"no shift" case were compared with those calculated where weighted averages 
were used in the time-shifting scheme. Table 3 displays the ratio of these 
two scatter functions by year for speeds (V) and densities (N). Ratios 
greater than 1 indicate that the scatter is reduced by time shifting using 
weighted averages. Also listed in Table 3 cure the yearly averages (as 
measured by ISEE 3) of the solar wind speed and density. We note that the 
importance of time shifting increases when solar wind speeds are low or when 
densities are high. Since higher solar wind speeds imply smaller t's, 
neglecting x (by not time shifting) would introduce greater error for 
low-speed intervals. These results also imply that variability exists in 
solar wind speeds even when the average speed is low. On the other hand, as 
evident by the density scatter function ratios, periods of high average 
density appear to be more variable than low average density intervals. 


TABLE 3. Scatter Function Ratios for ISEE 3 vs. IMP 8 Plasma Parameters 


a NO SHIFT ISE E 3 


Year 

# points 

a 

SHIFT 

(WEIGHTED) 

Averages 



V 

N 

V 

N 

1979 

1148 



390 

6.64 

1980 

1752 



378 

7. 24 

1981 

1542 

1. 16 

1. 10 

432 

6. 01 

1982 

1112 

1. 12 

2.00 

482 

8.68 
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Values of the actual plasma parameter functions are shown in Table 2. 
Density scatter functions determined using electron-based weighted averages 
are about 40% higher than those found when 5-min averaged ion-based densities 
were used* The velocity scatter functions differ by almost a factor of 3* 
These differences have three possible sources: (1) using electron-based 

plasma parameters in lieu of ion-based parameters, (2) constructing hourly 
averages by the method of weighted averages instead of collecting 5-min 
averages, and (3) the fact that the time periods compared are not the same. 
Regressing ISEE 3 hourly averaged ion-based densities and velocities to their 
IS EE 3 electron -based counterparts for 1978 through 1980 results in aZ of .17 
for log density and a Z of 33 km/sec for velocity based on 11,770 hours of 
overlapping data* Hence, we conclude that most of the differences between the 
density and velocity scatter functions can be attributed to the somewhat 
larger errors associated with electron-based parameters* 

In order to estimate what fraction of the irreducible differences results 
from the remoteness of ISEE 3, it is of interest to compare the scatter 
functions from Table 2 with those obtained in cross-normalizing pairs of 
near-Earth data sets. A special case is the pair of IMP 8 plasma data sets 
provided ty MIT and LANL. Virtually all of these data are from IMP 8 with the 
exception of some 1977 LANL IMP 7 data. Table 4 summarizes the scatter 
function data obtained upon cross -normalizing these data sets for 1977 through 
1984* A comparison of the Z functions to those of Table 2 suggests that 1/3 


to 1/2 

of the irreducible 

differences 

between 

ISEE 3 

and 

IMP 8 

V and 

log N 

values , 

and 

a yet smaller 

fraction for 

log T, 

may result 

from 

the remoteness 

of ISEE 

3. 

Rote the somewhat 

greater 

scatter 

in the 

1980-1982 

era. 




TABLE 4. 

MIT 

vs. LANL 

> Scattez 

■ Functions 

- IMP 

8 


Year 


Number of Hours 

Scatter Functions 


Z 

Functions 





V 

log N 

log T 


V 

log N 

log T 

1977 


4465 


4. 1 

.04 

.07 


5.8 

.05 

. 10 

1978 


2267 


4.9 

.04 

. 10 


6.9 

.06 

.14 

1979 


2213 


4.0 

.05 

.09 


5.7 

.06 

. 13 

1980 


1938 


4.5 

.06 

.10 


6.3 

.08 

. 14 

1981 


1575 


5.5 

.07 

.12 


7.8 

.10 

.16 

1982 


1235 


6.2 

.05 

.10 


8.7 

.08 

.14 

1983 


2400 


4.6 

.04 

.09 


6.6 

.06 

.13 

1984 


1616 


4.9 

.03 

.06 


7.0 

.04 

.11 
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In the original Interplanetary Medium Data Book of 1977, scatter 
functions for plasma and field parameters were listed for several pairs of 
data sets obtained from near-Earth spacecraft between 1963 and 1975. These 
scatter functions ranged between 7.5 and 17.0 km/s (V), 0.04 to 0.10 (Log-N), 
and 0. 06 to 0.20 (Log-T). There has been a trend towards decreasing random 
differences between data sets with time, due at least in part to the 
increasing number of instrument energy channels from whose count rates bulk 
flow parameters are determined. For magnetic field components, the scatter 
functions were typically in the 0.7 to 1.1 nT range, while for field 
magnitudes the scatter functions were either in the 0.3 to 0.6 nT range (8 
cases) or 0.9 to 1.1 nT range (4 cases). Unlike the IMP 8 plasma case, there 
has been no recent case of two magnetometers flown on the same spacecraft. We 

irreducibl « differences between the IMP 8 and corotationally 
s fted ISEE 3 hourly averaged IMF data sets are not significantly different 
statistically from the irreducible differences between non-time-shifted hourly 
data sets obtained for a pair of near-Earth spacecraft. In this sense, adding 
corotation-shifted ISEE 3 IMF data to our 1 AU, hourly average data 
compilation does not significantly diminish the reliability of the 
compilation. 


Cross -Normalization of Data Sets 


c u ^f 011814 ® 1 ^ next the systematic differences between individual data sets. 
Such differences are in contrast to tiie previously discussed random 
differences between data sets, and are reflected in the extent to which M best 
fit regression lines deviate from Y-X. Here the concern is to make the data 
sets contributing to this supplement as mutually consistent as possible, and 
separately to make these 1977-1985 data as consistent with earlier data as 
possible. Recall that while all parameters have been cross -compared in the 
past, only density and temperature have been normalized. 

Table 5 summarizes the ISEE 3 versus IMP 8 IMF regression results. In 
this Tabie, a is the scatter function. The "1-a range" column gives the range 
of P( ISEE3 ) over which the regression line lies within one sigma (parallel to 
ordinate) of Y-X. Typical field components and magnitude values are almost 
always deep within these ranges. Therefore, concluding that there would be no 
statistically significant gain in cross -normalizing the field data sets, we 
have not normalized them. 


TABLE 5. IMP 8 vb. ISEE 3 IMF Regression Parameters 


p 

a 

b 

P IMP " 

CL 

a + b Pisee 3 
1-cjt ranqe 


B x 

. 

o 

.98 

.97 

-66 to 70nT 

5-minute ion- based 

B y 

. 00 

.97 

1.01 

-47 to 47nT 

speeds used for x 

B z 

. 07 

1.03 

1. 12 

-55 to 50nT 

calculation 

B m 

-. 13 

1.00 

00 

V 

• 

O 

0 to nT 

1978-79 

B x 

.01 

.98 

1. 16 

-102 to 103nT 

hourly electron 

B y 

.02 

. 98 

1.27 

- 88 to 90 nT 

based speeds used 

B z 

• 02 

1.00 

1.33 

-388 to 380nT 

T calculation 

B m 

. 03 

.98 

.53 

0 to 40 nT 

1979-83 
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There are four plasma data sets to be cross-compared. Regression results 
for flow speeds, densities and temperatures are shown in Tables 6, 7 and 8 
respectively. In order to ascertain whether or not the regression parameters are 
time invariant, each comparison has been presented for each year. Also shown are 
the ranges of V, N and T over which the regression lines for V, log N and log 
lie within one sigma of perfect agreement (Y=X). 

Flow speed data generally agree to within a few km/sec. There are some 
yearly variations, but no trends are apparent. Even though for certain years the 
1-a range for the IMP(LANL) IMP (MIT) velocity regression is not ideal, we have 
elected not to normalize these data in order to avoid using a time varying 
normalization. A small but significant number of density and temperature values 
are measured in the parts of parameter space where the systematic differences are 
comparable to or greater than random differences. Thus, we shall follow our 
earlier approach of cross -normalizing the density and temperature data. 


ISEE 1 and ISEE 3 plasma data sets were compared to both the LANL and MIT 
data sets from IMP 8. Chaining of the derived regression equations demonstrates 
their mutual consistency. For example, the 1977 IMP (LANL ) vs. IMP(MIT) and the 
1977 IMP ( LANL ) vs. ISEE 1 relations of Tables 7 and 8 can be combined to yield: 

log N M i T = .16 + .84 log NiseE 1 

log T mit = .85 + .83 log Txs EE i 

These equations compare favorably to those found from the direct IMP (MIT) vs. 

ISEE 1 comparison. 

It remains to choose what normalization to apply to which data sets. Rather 
than presume to judge which of the four data sets is more likely to be absolutely 
correct, we shall normalize all densities and temperatures to IMP (LANL) values 
for historical consistency. In previous Data Books, the density normalization 
used for IMP (MIT) data from 1973 to 1978 was: 

log = .12 + .89 log N^it 

This equation is not statistically different from the normalization equation 
found by cross-comparing data from 1979 to 1984. Hence, we shall normalize the 
IMP ( MIT ) density data using the 1973 to 1978 relation. 

It is worth noting that although there are yearly variations in the 
normalization parameters, little is gained by applying a time varying 
normalization. The cr values are such that the 1-cr range calculated against the 
proposed normalization for any given year usually encompasses virtually all of 
the relevant parameter space. For 1980 and 1981 this is not true ^ since the 
IMP (LANL ) IMP (MIT ) regression line appears to undergo a "lowering". These 
"anomalous” years will be addressed subsequently in this section. 

The IMP (LANL) vs. IMP (MIT) temperature regression results for 1979 to 1984 
shown in Table 8 are essentially identical to the 1973-1978 results of 

log “ “’*62 + 1.11 log T^IT 

Hence for historical consistency the previous normalization will be used. Again, 
it is worth pointing out the year-to-year variations in the regression 
parameters. 
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TABLE 6. Flow Speed Regressions 


V LANL " a + 

b V MIT (for IMP 

8) 






Year 

# Points 

a 

b 

<3. 

1 - 

o range 

1977 

4465 

-7.0 

1.01 

4. 1 

100 

to 

1000 km/sec 

1978 

2267 

-1.6 

1.00 

4.9 

0 

to 

00 

1979 

2213 

-0.6 

0.99 

4. 0 

0 

to 

560 

1980 

1938 

0. 2 

0.99 

4. 5 

0 

to 

470 

1981 

1575 

4.3 

0.98 

5. 5 

0 

to 

470 

1982 

1235 

-2.7 

0.99 

6. 2 

0 

to 

1000 

1983 

2400 

1.2 

0.99 

4.6 

0 

to 

560 

1984 

1616 

1.6 

0.99 

4.9 

0 

to 

590 

1979-1984 

10,977 

0.3 

0.99 

4.9 

0 

to 

580 


V IMP - « + 

b V ISEE 1 








Year 

IMP 

# Points 

a 

b 

<L 

1 - a 

range 


Instrument 








1977 

LANL 

655 

3. 1 

1.01 

6.8 

0 

to 

670 

1978 

LANL 

635 

-2.0 

1.02 

11.7 

0 

to 

980 

1979 

LANL 

525 

-7.0 

1.03 

16.9 

0 

to 

1030 

1977-1979 

LANL 

1815 

-0.8 

1.02 

12. 1 

0 

to 

1060 

1977 

MIT 

461 

12.6 

0.99 

7.2 

220 

to 

2000 

1978 

MIT 

859 

12.8 

0.99 

9.5 

0 

to 

2300 

1979 

MIT 

921 

11.8 

1.02 

12. 3 

0 

to 

330 

1977-1979 

MIT 

2241 

17.0 

0.99 

11.0 

120 

to 

2700 


V IMP * a + 

b V ISEE3 







Year 

IMP 

# Points 

a 

b 

<5. 

1 - a 

range 


Instrument 







1978 

LANL 

815 

8.0 

0.99 

6.9 

0 to 

1200 

1979 

LANL* 

1148 

-27.2 

1. 07 

15.6 

80 to 

750 

1980 

LANL* 

1817 

21.4 

0.96 

24. 2 

0 to 

1300 

1981 

LANL* 

1605 

-27.3 

1. 08 

21. 1 

0 to 

710 

1982 

LANL* 

1177 

-21.8 

1. 06 

17.0 

0 to 

760 

1980-1982 

LANL* 

4599 

-6.3 

1. 03 

20.8 

0 to 

1190 

1978-1980 

MIT 

5008 

6. 0 

1.00 

8.4 

0 to 

oo 

1979 

MIT* 

1896 

-4.2 

1.03 

13. 1 

0 to 

810 

1980 

MIT* 

3097 

3.3 

1.02 

20.9 

0 to 

1100 

1981 

MIT* 

2794 

-19.5 

1.07 

15.8 

0 to 

580 

1982 

MIT* 

1930 

-17.7 

1.07 

19.6 

0 to 

670 

1979-1982 

MIT* 

9717 

-10.0 

1.05 

18.0 

0 to 

700 


♦Electron-based ISEE 3 parameters used* 
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TABLE 7. Density Regressions 


log n LANL = a + b( log % IT ) 


Year 

# Points 

a 

b 


1-a range 

1977 

4465 

.09 

.92 

.04 

2.7 to 41 

1978 

2267 

.07 

.94 

. 04 

1.3 to 150 

1979 

2213 

. 10 

.90 

. 05 

2.8 to 57 

1980 

1938 

.05 

.85 

.06 

0.6 to 7.9 

1981 

1575 

.08 

.82 

.07 

0.8 to 9.6 

1982 

1235 

. 18 

.85 

.05 

5. 1 to 51 

1983 

2400 

. 15 

.90 

.04 

10 to 145 

1984 

1616 

. 18 

.87 

.03 

12 to 55 

1979-198 

10,977 

.09 

.91 

.07 

0.9 to 110 


log = a + b(log NjggEl) 


Year 

IMP 

Instrument 

# Points 

a 

b 

Z. 

1-a 

range (N) 

1977 

LANL 

655 

.24 

.77 

.08 

3.6 

to 

30 

1978 

LANL 

635 

. 14 

.92 

.09 

1.7 

to 

1200 

1979 

LANL 

525 

. 16 

.88 

. 10 

1.5 

to 

320 

1977-1979 

LANL 

1815 

. 18 

.85 

.09 

2. 3 

to 

110 

1977 

MIT 

461 

. 19 

.84 

.06 

4 ‘§ 

to 

54 

14 

io iH 

1978 

MIT 

859 

.06 

.99 

.07 

10 4 

to 

1979 

MIT 

921 

.04 

.98 

.07 

10" 4 

to 

10® 

lx 

1977-1979 

MIT 

2241 

.07 

.96 

.07 

.2 

to 

10 4 


log N imp = 

a + b( log Njggg3 ) 

Year 

IMP 

# Points 


Instrument 


1978 

LANL 

815 

1979 

LANL* 

1148 

1980 

LANL* 

1817 

1981 

LANL* 

1605 

1982 

LANL* 

1177 

1980-1982 

LANL* 

4599 

1978-1980 

MIT 

5008 

1979 

MIT* 

1896 

1980 

MIT* 

3097 

1981 

MIT* 

2794 

1982 

MIT* 

1930 

1979-1982 

MIT* 

9717 


a 

b 

<L 

1-a 

range (N) 

19 

0.85 

.09 

3.6 

to 

136 

16 

0.88 

.07 

2.8 

to 

164 

07 

0.83 

. 10 

0.4 

to 

16 

05 

0.85 

.09 

0.3 

to 

15 

14 

0.82 

.07 

1.8 

to 

18 

06 

0.85 

.09 

0.4 

to 

20 

09 

0.93 

.07 

°‘? 

to 

740 

05 

0.99 

.06 

10 20 

to 

io oi 

00 

0.99 

.09 

io - 20 

to 

io 20 

01 

0.96 

.07 

.01 

to 

660 

09 

1.02 

.09 

.07 

to 

10 7 

00 

0.98 

.08 

10“ 5 

to 

10 5 


♦Electron-based IS EE 3 parameters used. 
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TABLE 8. Temperature Regressions 


log t LANL = a + b( log T MIT ) 


Year 

# Points 

a 

b 


1-a range 

(T) 

1977 

4465 

-.77 

1. 15 

. 07 

2.0x10 4 

to 

4. 6x10 ! 

1978 

2267 

-.75 

1. 16 

. 10 

8. 3x10 3 

to 

6. 6x10 

1979 

2213 

r- 

00 

. 

i 

1. 18 

.09 

1.5x10 4 

to 

3.9x10 

1980 

1938 

-1. 14 

1. 24 

. 10 

1. 3x10 4 

to 

1.9x10 

1981 

1575 

- .43 

1. 11 

.11 

3. 5x10 2 

to 

3. 7x10 

1982 

1235 

.25 

0.96 

. 10 

7. 1x10 2 

to 

4. 5x10 

1983 

2400 

- .72 

1. 15 

.09 

I.OxlO 4 

to 

5. 3x10 

1984 

1616 

- .04 

1.02 

. 08 

5.0x10" b 

to 

1. 3x10 

>79-1984 

10,977 

- .62 

1. 13 

. 10 

4. 4x10 3 

to 

5. 2x10 


°K 


log T imp = a + b( log T ISEE i) 


Year 

IMP 

Instrument 

# Points 

a 

1977 

LANL 

655 

.21 

1978 

LANL 

635 

.30 

1979 

LANL 

525 

-.97 

1977-1979 

LANL 

1815 

. 14 

1977 

MIT 

461 

.81 

1978 

MIT 

859 

.80 

1979 

MIT 

921 

.46 

1977-1979 

MIT 

2241 

.85 


b 

c x. 

1-a range 

(T) 

0.95 

.06 

1 . 7x10 2 

to 

2. 1x10 

0.93 

.09 

2. OxlO 2 

to 

8.7x10 

1. 17 

. 14 

3. 6x10 4 

to 

1. 2x10 

0.96 

.11 

M 

to 

3. 4x10 

.84 

.08 

1 . 8x10 ^ 

to 

4.9x10 

.83 

. 11 

6. 3x10* 

to 

3.7x10 

.86 

. 14 

10 3 

to 

4.6x10 

.80 

. 14 

2.5x10 

to 

2. 2x10 


log t IMP = 

a + b( log Tisee3) 




Year 

# Points 

a 

b 


1978-1980 

1978 

5008 (MIT) 
815 (LANL) 

-.06 

.34 

1. 04 
0.94 

. 08 
.06 


1-a range (T) 

.02 to 1.6x10 5 
3. 3x10** to 2. 2x10 7 
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Cross-comparisons of ISEE 1 with either IMP instrument were performed on 
a fairly small number of points. Thus the significance of the year-to-year 
variations is not easily determined. Chaining the IMP(LANL ) -IMP (MIT ) with the 
IMP (MIT) -ISEE 1 equations yields results consistent with the direct IMP (LANL) 
ISEE 1 equations for both densities and temperatures. Temperatures will not 
be normalized since the 1977-1979 regression parameters do not differ 
significantly from Y=X. We do normalize ISEE 1 density as indicated in Table 9. 

A cross -comparison of ion-based ISEE 3 densities with IMP (MIT) densities 
was done for 1978 through day 48, 1980 (Table 7). Similar regression 
parameters are obtained when ISEE 3 electron-based densities are used for 


The IMP(LANL) ISEE 3 density regression parameters for 1979 and 1982 lie 
within a 1-a range of the normalization used in Supplement 2 (log N^^ = .20 
+ .83 log N isee 3 ). However, 1980 and 1981 IMP(LANL) ISEE 3 regression 
parameters have a different character. The regression line appears to be 
lowered in a similar manner to that found from the IMP(LANL )-IMP(MIT ) 
cross-comparison. Since IMP (MIT) ISEE 3 regression parameters do not exhibit 
any anomalies for 1980 and 1981, we assume that any instrumental time 
variability is in the IMP (LANL ) instrument. Thus we normalize all IMP (MIT) 
and ISEE 3 densities using the same relation that was used in Supplement 2. 

As explained earlier, IMP(MIT) parameter values are selected preferentially 
over IMP ( LANL ) values when both are available for any given hour. Because of 

this, there is only a small number (~300) of 1980-1981 hours having IMP (LANL) 
data. 


To further support the ignoring of time variations in the normalizations, 
we depict in Figure 12 the yearly averages of the densities, velocities, and 
temperatures for the two IMP plasma instruments. These averages are based on 
the overlapping hours listed in the first part of Tables 6, 7 and 8. 

Velocities agree quite closely. Average temperatures are not as close, but 
LANL temperatures are consistently higher than MIT temperatures. The average 
LANL densities are lower than MIT densities for 1980 and 1981 yet higher for 
all other years. 

Table 9 summarizes the normalizations used in this book. 


TABLE 9. Normalization Parameters for N and T 


P NORM = a + b Pps 


P = log N 

JDS a b 


P = log T 
a 


b 


IMP-LANL 

0 

IMP-MIT 

. 12 

ISEE 3 

. 20 

ISEE 1 

. 18 


1.00 

0 

1.00 

0.89 

-.62 

1. 11 

0.83 

-.55 

1.07 

0.85 

0 

1.00 
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500 

475 


y 450 
(km/sec) 425 
400 
375 


(10 4o K) 
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Data Coverage 


The percent of coverage of the composite data set over the 1963 to 1985 
time period is shown in Figure 13. Of the 188,568 hours included in Bartels' 
solar rotations 1783 to 2073, there are 98,771 hours with field and plasma 
data of which 72,548 hours have field and plasma data from a common 
spacecraft, 22,031 hours with field data only, 24,960 hours with plasma data 
only, and 42,806 hours with no interplanetary plasma or field data. 

Data Presentation 


This third Data Book Supplement consists of graphical (Volume 3) and 
tabular (Volume 3A) presentations of some of the parameters of the composite 
data set. In Volume 3, there cure two plots for each solar rotation in which 
any plasma or field data were obtained. On facing pages , for convenience in 
lining up features in the data, are found a plot of plasma data (bulk speed, 
density, and proton temperature) and a plot of field data (average magnitude, 
geocentric solar magnetospheric (GSM) B z component, and geocentric solar 
ecliptic (GSE ) latitude and longitude angles of the average field vector). 

Note that on those rare occasions when the parameter values exceed the allowed 
range, a heavy mark is placed near the edge of the plot. For such cases the 
reader is advised to consult the data listings (Volume 3A) for appropriate 
numerical values. 


Additional Data Availability 

In addition to the parameters listed and plotted herein, the data set 
from which this Data Book Supplement was generated also contains additional 
IMF parameters (e.g. , By, Bz in GSE coordinates), additional plasma parameters 
(flow direction), standard deviations in IMF and plasma parameters, and 
geophysical and solar activity indices (Kp, C9, Dst, R). 

This data set is available both online, on the NSSDC VAX, and on magnetic 
tape. The data set is updated as NSSDC receives additional relevant data, at 
a typical frequency of every several months. New hardcopy Supplements are 
issued only every few years. 

A word on day-numbering conventions is appropriate. When the first 
OMNItape was created, most input data used the convention that January 1 is 
Day 0. This convention was employed for the OMNItape, and has been continued 
for all subsequent tape versions. However, it is recognized that this is a 
minority convention. Accordingly, the Data Books and the recently created 
online "OMNIfile" both use the convention that January 1 is Day 1. 

In all versions of this data set, missing parameter values sure filled 
with zeros. 
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The Online File 


The data set is online from 1976 onward. Earlier portions can be brought 
online in response to demand. The data set may be accessed over SPAN by $SET 
HOST NSSDC , USERNAME =NSSDC • This interface gives a menu of NSSDC online 
services, one of which is access to the "OMNIfile." Currently, the user may 
view the file format and may select and list at his terminal any subset of 
parameters for any days of interest. Additional capabilities, such as linking 
to an NSSDC -supplied READ subroutine to read OMNIfile records, and creating a 
subfile for downloading to the user's VAX, are currently being developed. 

The Magnetic Tape 

ASCII or IBM/binary magnetic tape versions of the data set from 1963 
onward are also available. Copies of either of these tapes (or a reformatted 
version), with documentation, may be ordered electronically from another menu 
option of the USERNAME =NSS DC account discussed in the previous paragraph, or 
by request to: 

National Space Science Data Center 
Code 633.4 

NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 
Telephone: (301) 344-6695 

Telex No. : 89675 NASCOM GBLT 

TWX No. : 7108289716 

Researchers who reside outside the United States would contact: 

World Data Center A for Rockets and Satellites 
Code 630.2 

NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 U.S.A. 

Telephone: (301) 344-6695 

Telex No. : 89675 NASCOM GBLT 

TWX No. : 7108289716 
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